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O , ABSTRACT 

' Destriping methods for constructing maps of the Cosmic Microwave Background 

(CMB) anisotropies have been investigated extensively in the Uterature. However, 
their error properties have been studied in less detail. Here we present an analysis 
of the effects of destriping errors on CMB power spectrum estimates for /'/ancfe-like 
' scanning strategies. Analytic formulae are derived for certain simple scanning geome- 

tries that can be rescaled to account for different detector noise. Assuming Planck-Vike 
low-frequency noise, the noise power spectrum is accurately white at high multipoles 
{£ }t 50). Destriping errors, though dominant at lower multipoles, are small in com- 
parison to the cosmic variance. These results show that simple destriping map-making 
, methods should be perfectly adequate for the analysis of Planck data and support the 

' arguments given in an earlier paper in favour of applying a fast hybrid power spectrum 

, estimator to CMB data with realistic 'l/f noise. 

^ . 
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INTRODUCTION 



The problem of constructing a map of the CMB anisotropies from a set of time-ordered data (TOD) has been studied by 
_ ' many authors. The methods can be broadly divided into two classes: 'optimal' methods that provide a least squares map 
^ ! making solution {e.g. Wright 1996; Wright et al. 1996; Tegmark 1997a,b; Borrill et al. 2001; Natoli et al. 2001; Dore et 
$_( 1^- 2001) and approximate 'destriping' methods {e.g. Burigana et al. 1997; Delabrouille 1998; Maino et al. 1999; Revenu 



et al. 2000; Keihanen et al. 2003). A brute force application of 'optimal' methods requires the inversion of large matrices 
and is computationally impractical for large TODs such as those expected from WMAP and Planck* . As a result, iterative 
algorithms have been developed {e.g. Wright et al. 1996; Natoli et al. 2001; Dore et al. 2001) which do not require matrix 
inversions. Nevertheless, for P/ancfc-sized datasets, these iterative algorithms are computationally expensive and require the 
use of supercomputers. 

Destriping algorithms are well suited to a Planck-type scanning strategy in which the sky is scanned many times on rings. 
The TOD can then be averaged on rings and the effects of low frequency noise noise approximated by a constant offset for each 
ring. The offsets can be determined from the ring overlaps. Destriping algorithms are conceptually simple and computationally 
fast. Even for Planck-sized TODs it is practical to apply destriping map making methods on many thousands of simulations 
to test the effects of various systematic errors (see e.g. Poutanen et al. 2004). 

The motivation for this paper is twofold. Firstly, although a number of authors have investigated destriping algorithms, 
almost all of this work has been numerical. A notable exception is the paper by Stompor and White (2004) who present an 
analysis of destriping errors for some simple scanning strategies. One of the aims of this paper is to develop on the work 
of Stompor and White and to derive an analytic model of the effects of destriping errors on the CMB power spectrum for 
Planck-like scanning strategies. This is useful because it helps in developing an understanding of the map-making process 
and how the errors depend on the parameters of the experiment. This analysis also sheds light on the differences between 



* See the Planck web-site http://www.rssd. esa.int/index.php?project=PLANCK. 
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destriping and 'optimal' map-making methods. In particular, whether the extra computational cost and complexity of an 
optimal method actually produces any significant improvement on simple destriping. 

The second motivation for this paper follows from the need to compute estimates of the CMB power spectrum, C't, 
rapidly and acurately. In an earlier paper (Efstathiou 2004, hereafter E04) a fast hybrid estimator was developed that 
combines a quadratic maximum likelihood estimator at low multipoles with a set of 'pseudo-C^' estimates at high multipoles 
with different pixel weightings. In E04 this method was tested against numerical simulations that used a realistic scanning 
strategy for Planck, but assumed uncorrelated white noise. In this approximation, the hybrid estimator was shown to be 
very close to optimal and, importantly, can provide an accurate estimate of the covariance matrix (ACeACei). However, in a 
realistic experiment, striping errors will introduce correlations in the noise. The question then arises as to whether the hybrid 
estimator is applicable, for example, is there a natural angular scale below which the noise can be assumed white, and if so, 
what fixes this scale? 

Our goal, therefore, is to investigate the effects of destriping errors on the CMB power spectrum for realistic scanning 
strategies and noise models. It should be emphasised that we do not attempt to develop a new map making technique, nor 
to investigate 'real world' complexities such as asymmetric beams, positional errors or non-stationary noise, though we will 
comment on how some of these aspects may effect our results. 



2 OVERVIEW OF MAP MAKING WITH DESTRIPING 

We denote the noise contribution to the TOD by n{t), which is considered to be a vector specified at integer values of the 
sampling frequency tsamp, i-e. Hi = n(ifsamp), i = 1, ■ ■ ■ N. The discrete Fourier transform of the noise TOD is denoted h{i') 
and the power spectrum (|n(i')|^) is assumed to be of the form 

P„(^) = (i + V i<^.<^_, (1) 

r T - - 2t 

samp 

where i^knee is the 'knee frequency'and T is the total length of the time-stream. The variance of the noise contribution given 
the power spectrum (1) is 

{\n{t)f) =alil + !^k„eeln(0.5r/i,amp)) , (2) 

and thus cr^ fixes the overall amplitude of the noise and is equal to the variance of white noise in the limit i^knee — > 0. 
The actual TOD is the sum of the true sky signal Sp and the noise TOD 

where Pip is a pointing matrix mapping a pixel p on the sky to the pixel i in the TOD. For a Planck-typc scanning strategy, 
the same circle on the sky is mapped A^ropoint = 60 times before the satellite is repointed. We will denote the spin period by 
tspin and the repointing time interval by trepoint (= Nrepointtspin) ■ As explained in the introduction, the main aim of this paper 
is not to develop an 'optimal' map-making technique, but rather to gain an intuitive understanding of the error properties of 
simple destriping map making methods. We therefore assume perfect pointing and simply average the TOD on each scanning 
ring: 

~ *i-Hj^*fringpix-Kfc-l)JVrepomt''ringpix' .^ringpix = tspin/tsamp- (4) 

3=0 

where ATrmgpix is the number of pixels withing a single ring pixel and Xi is the mean vaue of the TOD in ring pixel i of ring k. 

This simple averaging is the maximum likelihood solution for map making on a ring if the pointing is assumed to be perfect. 
It is possible to improve on (4) to account for imperfect pointing by solving the maximum likelihood equations to reconstruct 
a map from all scanning rings within a single pointing period (see, for example, van Leeuwen et al. 2002). 

For the simulations presented here, we adopt the parameters listed in Table 1 unless stated otherwise. Table 1 also serves 

as a summary of the notation used in this paper. The knee frequency in Table 1 has been chosen to be representative of the 
70 GHz channel of the Planck Low Frequency Instrument (LFI) (see Tuovinen 2003) . This is an interesting case because 
the knee frequency is about twice the spin frequency. The knee firequencies for the Planck High Frequency Instrument (HFI) 
should be smaller than the spin frequency. This case is less interesting because if t'knoo ^ i^spin, it is a very good approximation 
to model the low frequency noise as a constant offset in each ring. The noise level dn for the simulations has been chosen so 
that the CMB power spectrum is noise dominated at multipoles t ^ 300 (rather than to match the noise for any of the Planck 
detectors). The input CMB power spectrum, Ci, is that of the concordance ACDM model favoured by WMAP (Spergel et 
al. 2003). As in E04, unless stated otherwise beam functions will not be written explicitly in equations, thus Ce will ususally 
mean Cibl, where hi is the spherical transform of a symmetric Gaussian beam. The remaining parameters, such as the number 
of rings, number of ring pixels, map pixel size etc were chosen so that large numbers of simulations could be run quickly. 
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Table 1: Notation and parameters 



symbol description value 





number of rings 




2160 


^^ringpix 


number of pixels per 


ring 


2160 


rspoirit 


number of rings per repointing 


60 


T 


length of TOD 




7.78 Ms 




pointing period 




3600 s 




spin period 




60 s 


^samp 


sampling period 




27.78 ms 




spin frequency 




16.67 mHz 


^knee 


knee frequency 




30.00 mHz 


t^max 


maximum frequency 




18.00 Hz 


(TiV 


noise amplitude 




1.36 mK 


■^^ring 


ring width 




10 arcmin 


^^ringpix 


size of ring pixel 




10 arcmin 


^^map 


size of map pixel 




15 arcmin 


^beam 


beam width FWHM 




30 arcmin 


-^map 


number of map pixel: 


3 


659676 




signal+noise in TOD 


pixel i 






mean signal+noise in 


pixel i of ring k 






offset of ring k 






trip 


map pixel p 







The basic assumption behind destriping techniques is that low frequency drifts in the TOD can be accounted for by 
adding a constant offset tk to each ring. The ring offsets can be determined by minimising 

5= 5^(£fc+5*-q-S^-f + A j^^q J , (5) 

where the notation j <Z i indicates that the ring pixels j and i overlap the same map pixel nip. The second term in equation 
(5) is included to enforce the condition = 0. The offsets are therefore given by the solution of the linear equations 



Y^{tk + x'l -ti-x])Jt-\\y^ei\ =0. (6) 

ijCi \ I J 

There has been some discussion in the literature concerning the weighting of the term in the first summation in equation (5). 
Equation (5) assigns equal weight to each overlapping pixel, as in Maino et aZ.(1999). Delabrouille (1998) assigns a weight 
w = l/{np — 1), where rip is the total hit count in map pixel p, while Keihanen et al. (2003) assign a weight w = l/up. 
For Planck-Vike scanning strategies the differences between destriping using these weight functions are much smaller than the 
striping errors themselves (see Figure 3 of Keihanen et al. 2003) and so the weight function will be set to unity throughout 
this paper. 

Equations (6) can be solved by a matrix inversion and the solution is independent of the regularising parameter A, 
provided it is chosen to be large enough. Alternatively, these equations can be solved by iteration by setting 

^k»^Y.^x]-x\), (7) 

IjCi 

adding the derived offsets to Xi and re-evaluating equation (7) until the offsets (7) converge to zero. In (7), Nok is the number 
of pixels on all rings that overlap pixels on ring k. This algorithm is identical to the plate matching procedures applied to 
create galaxy catalogues from photographic data (Groth and Peebles 1986; Maddox, Efstathiou and Sutherland 1996). If the 
number of overlaps per ring is large, then an accurate estimate of the variance of the offsets can be derived from the first 
iteration of equation (7), 

~ i\t2 at ^ AT ' ^ ' 

ok ^^repomt ok -^^repoint 

where NoWck is the number of identical ring pixels summed over all pairs of rings I, I' that overlap with ring k. In equation 
(7) we have assumed that the noise in each averaged ring pixel is white, with variance crJ;/A'repoint, which is accurate on a 
single ring even in the presence of 1// noise unless i^nee ^ fspin. The second term on the right hand side of equation (8) 
applies for a Planck-type scanning strategy for which NoWck ^ ^ofc- 
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Figure 1. The correlation function of the ring offsets for two Planck-type scanning strategies; (a) for a bore-sight angle di, = 85° and 
spin axis aligned with the ecliptic plane; (b) for Of, = 85° and with a slow sinusoidal precession of = 5° above and below the ecliptic 
as discussed in the text. The indices i and j refer to the ring number. 



We consider simple Planck-type scanning strategies with the spin axis either ahgned with the ecliptic plane, or with a 
slow precession of 5° sin(2c/>e) about the ecliptic plane, where (t>e is the ecliptic longitude. The sky is scanned with a single 
detector at a' bore-sight' angle of 9i, with respect to the spin cixis. After a complete uniform sweep of the ecliptic plane, the 
time-stream maps to a set of A'rings each of angular width Airing. 

Figure 1 shows the correlation functions of the ring offsets averaged over 250 simulations for each of the two scanning 
strategies discussed in the previous paragraphs. There figures indicate the following: 

(a) The dispersion in the ring offsets is {eiy/^ « 11 nK^ in excellent agreement with equation (8). This is much smaller 
than the white noise level of 176 /uK on a single ring because the number of overlaps for a PZarecfc-type scanning strategy is 
large. As emphasized by Stompor and White (2004) the pixel noise of the resulting maps will be predominantly white and 
uncorrelated. 

(b) It is interesting to compare the ring variance with the expected signal variance for the case of no precession: 

{el) = ^ Y,^2£ + l)Ce [Pe{cos e^)f = (14.2 i^Kf. (9) 

Thus, for the parameters adopted in this paper, the oflfset variance arising from ring pixel noise is comparable to the true 

signal variance. 

(c) The offset correlation functions in Figure (1) drop rapidly to zero. To high accuracy, we can model the striping errors 
as a set of random offsets with dispersion a and effective ring- width Aa as illustrated in Figure 2. As shown in the next 
Section, The contributions of these errors to the spherical harmonics and power spectrum of the destriped maps can, in this 
approximation, be calculated exactly for the case of a perfect ring torus. 

(d) There is almost no perceptible difference in the correlation functions for the two scanning strategies. We would therefore 
expect (and this is verified in Section 4) that the effects of destriping errors on the power spectra should be very similar for 

these two scanning strategies. 

(e) The striping errors arising from the dispersion in the ring offsets are 'irreducible' errors. By this, we mean that these 
errors are fixed by the instrumental white noise on the rings and the crossing points (interconnectedness) of the TOD. They 
cannot be reduced by applying more time-consuming 'optimal' map-making methods (see Section 5). However, since the knee 
frequency i^knee exceeds both the repointing frequency 1/irepoint and the spin frequency, the averaged ring data Xi will contain 
low amplitude gradients associated with '1//' noise. These gradients can, in principle, be removed by modifying the destriping 

t For Planck the dispersion will be smaller because the detector noise is smaller than assumed here 
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code to determine from the crossing points a few low order coefficients in, say, a Legendre poynomial or Fourier expansion, 
(Delabrouille f 998; Keihanen et al. 2003) However, as we will show in Section 3, the effects of these gradients are smaller than 
the effects of the offset errors. 



3 ANALYTIC MODELS OF DESTRIPING ERRORS 

In this Section we consider a scanning strategy which leads to a perfect ring torus, i. e. the spin axis is aligned with the ecliptic 
plane as the sky is scanned by a single detector with a bore-sight angle of 9i,. For such a scanning strategy, the distribution 
of hit-counts on the sky will follow a distribution with ecliptic latitude of 

r 7V(1 -cos2 6l/sin^6'6)~^/^sin6'd6l, cos < sin 06, 

{ 0, cose> sin 06. ^ ^ 

The hit-count distribution is therefore highest at cos6' — sin 6^6 and lowest at the ecliptic 8 = tt/2. 

As described in the previous Section, the offsets et of each ring can be modelled as a set of independent Gaussian random 
variates with dispersion a. The map constructed from these ring offsets (which we will refer to as the error map) will contain 
most of the information on pixel correlations introduced by the map-making process. Our goal in this Section is to compute 
the effects of these errors on the CMB power spectrum. 

The spherical harmonic transform of the error map can be written as 

dim = WikATikilikYlrn{9ik, (t)ik), (11) 

ik 

where the index k denotes the ring number, i denotes the pixel number within the ring and Qik is the solid angle of the ring 
pixel. The weight factors Wik account for the averaging of the ring pixels in constructing the map and thus are proportional 
to the inverse of the hit count distribution of equation (10). 

To evaluate equation (11), reorient each ring to a new coordinate system {6' , (j)') in which the spin axis is aligned with 
the new z' axis. The spherical harmonic transform for each ring is then, 

/efc 1 mm Z"^" / im ' 

a'lm = 7: <^kA'^ Pji^ (cos 6b) sin 9b Aa I | sin (/i'le""*^ d<jf>', (12a) 
2 Jo 

where ak is the constant ring-offset ATik = Uk and the A™ are the normalising factors of the spherical harmonics 
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-[ 47r {£ + my.) ' 
Notice that the hit count distribution has been normalised so that the sum of the weight factors 

Wikflik = ^TTsmOb (13) 

ik 

for a complete ring torus (which completely covers the sky twice for the case di, = n/2). Performing the integral over 0' in 
equation (12a), 

a'em = cTk sin ObAaATPPicos 6b) J- — - — ^(1 + (-1)™) = crksm6bAaK^{coaBb), (14) 
where K'^ is zero for all odd values of m. Transforming back into the original coordinate system 

o|m = ^ ^ a't^iDini {ak,Pk,7k), (15) 

k m' 

where the D^^^ are the Wigner D-matrices (see e.g. Brink and Satchler 1993; Varshalovich, Moskalev and Khersonskii 1988) 
and a, (3 and 7 arc the Euler angles relating the two coordinate systems. In our case, Uk = 0, /3fc = 7i"/2, hence in terms of 
the real reduced rotation matrices (15) is 

o|m = ^^crfcAasin6'6-KT d^/^(7r/2)e"'™''"= . (16) 

k m' 

Defining the power spectrum of the error map as 
C'l = ^ 1) 53(«?ma?m, ) (17) 

m 

(where the tilde on C signifies that the power spectrum is computed from the af^ coefficients computed on the incomplete 

sky if Ob 7^ 7r/2) and using the relation 

^dm'm(T/2)dm"m(7''/2) = Sm'm", (18) 
m 

we find 

C! = ^^l^c^'Aasin^ 9b J2 l^r(cos^b)|^ (19) 



For the special case 9b = 7r/2, the factor \Kp\ is proportional to Pp{0)^ and so vanishes for odd values of {£ — m). Since 
is zero for odd values of m, it follows that (7| = for odd values of I. An alternative derivation of the C| for the special 

case 9b — n/2, (in which the spherical harmonic transform for each ring is evaluated in a coordinate system with the ^;-axis 

perpendicular to the spin axis) gives 

C*! = ^a^Aah, (20a) 
where 

f - r6icos6'P^(cos6')d6' £ even, e^O 
^^ = \ 0, ^odd. 

For large values of £ the integral in equation (20b) can be approximated by 

poo 

9cos9Pe{cosO)dO « tt / Jo((^ + 1/2)0) = , (21) 



Jo 



'0 

hence 



; + i/2)' 



C! « j^^a^Aa, £ even. (22) 

Equation (22) is, in fact, an extremely good approximation (to within 2%) to the exact answers of equations (19) and (20a) 
for values of £ as small as £ = 4. 

Some examples of the error power spectra for three values of 6b are shown in Figure 3. For each value of 6b, 10* simulations 
were generated each with 720 rings assigned random offsets. The rings were then mapped on to the igloo pixelization scheme 
described in E04 with 0.25° x 0.25° pixels simply by assigning the nearest map pixel to each ring pixel and averaging over all 
of the ring pixels assigned to any particular map pixel. The power spectra were computed from the igloo maps using using 
fast spherical transforms. As can be seen from Figure 2, the finite sizes of the ring and map pixels introduce some structure 
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Figure 3. The power spectrum of destriping errors for ring tori with various values of 6;,. In each figure the filled (blue) points show the 
results from numerical simulations and the solid (red) lines show the analytic expression of equation (19). The dot-dashed (green) line 
shows the simple analytic approximation of equation (22). 
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in the final error maps. Nevertheless, the mean power spectra for the error maps (shown by the points in Figure 3) agree 
perfectly with the analytic results of equation (19) which were derived in the continuum limit. Evidently, the effects of finite 
pixolisation and ring widths arc negligible and hence the analytic model developed in this Section gives an extremely accurate 
representation of the destriping errors. 

If the primordial fluctuations are Gaussian, the spherical harmonic coefficients will satisfy 

(«lmfl|*m') = CiSu'Smm'- (23) 

The striping erros will, however, introduce correlations in the acm- Prom equation (16) it is straightforward to show that the 
correlations introduced by striping arc given by 

{at^atr^,)=27Ta\m''eiAa6^m' ^ Kp Kpdl^^{7T/2)di^^,i7r/2). (24) 
For the special case 0b = n/2, these correlations can be written as 

{aemapm'} = 27rcT^Aa S^rim'^rE^ (1 + (-1)™), (25a) 
where 

Er= [ ATPrilAdlJi. (25b) 



If the vector is ordered as (m.l), {i.e. (0,i?min), ■ • • (0,i'max), (l,^min), • • • (1,-^max), etc) the covariancc matrix {(iim(^t*m') 
will have a block diagonal structure. Figure 4 shows examples of these covariance matrices for the three scanning strategies 
discussed above. As we will show in the next Section, these ofF-diagonal correlations will be much smaller than the diagonal 
components for the parameters of a realistic experiment. However, since the departures from Gaussianity in most inflationary 
models are expected to be small (for a comprehensive review see Bartolo et al. 2004), residual map-making errors could prove 
problematic in testing for non-Gaussianity. For the remainder of this paper, we will concentrate on errors on estimates of the 
CMB power spectrum, though it is important to recognise that oven if map-making errors can be shown to have a negligible 
effect on the power spectrum, they may be important for other statistical tests. 



4 SIMULATIONS WITH REALISTIC NOISE 

In this Section, we describe the results from a set of simulations with realistic '1//' noise (equation 1) with the parameters 
given in Table 1. As explained in Section 2, the resolution and sizes of the maps and ring sets were chosen so that large numbers 
of simulations could be run quickly while demonstrating the salient features of the map-making problem. As a further speed 
up, 1//- noise was generated by using an FFT for frequencies below 0.133 Hz. Above this frequency the noise was assumed 
to be white, which has the additional advantage that the white noise can be added to ring pixels 'on the fly' so that it is 
never necessary to store a complete TOD in memory. A complete simulation, including noise generation, destriping and power 
spectrum estimation takes approximately 70 seconds on an single 1.4 GHz Itanium 2 processor. Splitting the noise into a low 
frequency '1//' component and a high frequency white noise component has the additional advantage that the effects of low 
and high frequency noise on destriping can be investigated separately. 

We have run two sets each of 250 simulations for the two scanning strategies adopted for Figure 1, namely, Ob = 85° 
with no precession and 6b = 85° with a slow sinusoidal precession of ±5° . In each case, we generated three ring-sets for each 
simulation: one with white noise only, one with low frequency '1//' noise only, and one using the sum of those two noise 
models. These ring-sets were passed through the destriping algorithm to produce three maps per simulation. The averaged 
power spectra for the two sets of simulations are plotted in Figure 5. 

The upper panels in Figure 5 show the power spectra for the case of white noise only on the rings. This is the case 
that is closest to the analytic model of the previous Section. As expected, the analytic model of equation (19) summed with 
the appropriate constant white noise level, provides an excellent match to the simulations when proper allowance is made to 
calculate the effective ring width using the correlation functions of Figure 1. This is true even for the slow precession scanning 
strategy shown in Figure 5d. The main effect of a the slow precession is to fill in the coverage gaps at the eclipic poles, but 
the ring pattern is so close to a perfect ring torus that the residual striping errors produce a nearly identical effect on the 
power spectrum. 

Figures 5b and 5e show the residual striping errors when only low frequency '1//' noise is included. These errors arise 
primarily from residual gradients on the rings and so their amplitude depends on the knee frequency. Had we adopted a knee 
frequency much smaller than the spin frequency, these errors would have had a much lower amplitude. Nevertheless, even for 
the parameters adopted here, the amplitude of these errors is considerable smaller than the errors caused by the dispersions 
in the ring offsets. Notice that these errors also decay roughly as l/£, as for the error power spectra for pure white noise. 
As mentioned earlier, Delabrouille (1998) and Keihanen et al. (2003) have explored fitting low order functions for each ring 
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h h 

Figure 4. The covariance matrices for foi" tlie simulations shown in Figure (3) compared to the analytic expression of 

equation (24). The indices ii and 12 are ordered as (m, i), {i.e. (0, 0), . . . (0, ^max), (1, 0), . . . (1, ^max), etc) with m and £ running from 
to 10. 
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Figure 5. Averages over 250 simulations of noise for a single detector with 9i, = 85° and no precession (Figures 5a, b, & c) and with a 
slow sinusoidal precession of 5° amplitude as discussed in the text (Figures 5d, e & f). Filled circles in the upper panels (a &: d) show 
residual errors after destriping assuming only uncorrelated white noise on the rings. Filled circles in the middle panels (b & e) show 
residual errors after destriping assuming only low frequency 1/f noise. Filled circles in the bottom panels (c & f) show residual errors 
including white noise and low frequency 1/f noise. The solid lines in the upper panels show the model for destriping errors plotted in 
Figure 3b together with the white noise level appropriate to the maps. The solid lines in the lower panels show the same destriping model 
renormalised to fit the simulations. 
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during destriping, rather than a single offset efe. The results are mixed and in fact Keihanen et al. (2003) find that including 
more parameters actually produces larger power spectrum errors than simply fitting constant offsets. This is not suprising 

because the errors on the individual crossing points arc dominated by the white noise. If the white noise level is high, and the 
knee frequency is low, there may not be enough crossing points to determine more than a single constant offset per ring with 
any precision. 

The points in the lower panels in Figure 5 show the errors for the full noise model. These are just the sum of the white 
noise and '1//' errors. The theoretical models, plotted as the solid lines, are simply equation (19) rescaled to provide a good 
match to the simulations, summed with the appropriate constant white noise level. These provide an excellent match to the 
simulation results. The general shape of the power spectrum errors C| plotted in Figures 5c and 5f is well known from previous 
numerical work on destriping (Delabrouille 1998; Maino et al. 1999; et al. 2000; Keihanen et al. 2003). However, the results 
presented here explain why the errors have this particular form and how they depend on the parameters of the experiment. 

Finally, in Figure 6 we illustrate the effects of '1//' noise on a simulated map of the CMB sky for the scanning strategy 
with a slow sinusoidal precession. Figure 6a shows the power spectra of the input map, the destriped map and for the difference 
map. The power spectrum of the difference map is shown on a greatly expanded scale in Figure 6b, together with the analytic 
model plotted in Figure 5f. Figure 6c shows the differences between the power spectra of the input and output maps but with 
the white noise level subtracted, again plotted on an expanded scale. The solid lines show the expected errors 



using the error model of Figure 5f. The dotted lines in this Figure show the expected error from white noise alone. Notice that 
pure white noise is an excellent approximation to the errors for i ^ 10 and that even at low multipoles the destriping errors 
are much smaller than the cosmic variance (AC|) — 2Cf /(2£+ 1). For example, the error in the quadrupole amplitude for this 
realisation is AT2 = £{£ + l)AC2/27r = 10.4 /iK^, compared to the expected cosmic variance of ATf = 717.8 nK^. Compared 
to the cosmic variance, the destriping errors at low multipoles can be ignored. Furthermore, errors of this magnitude are much 
smaller than the errors of 50-100 fxK^ expected from inaccurate subtraction of the Galaxy (c/the discussion of the effects of 
Galactic subtraction on the low multipoles measured by WMAP, Bennet et al. 2003; Slosar and Seljak 2004). 

5 CONCLUSIONS 

In this paper, we have presented an analytic analysis of the effects of destriping errors on the CMB power spectrum for 
various scanning strategies. Destriping errors produce a characteristic error power spectrum of the form shown in Figure 
3 which decays roughly as l/£. The amplitude of the error power spectrum is determined by the number of intersections 
in the ring set which depends on the scanning strategy. In agreement with Stompor and White (2003), and with earlier 
numerical work, there are so many interconnections in a Planck-like scanning strategy that the destriping errors should be 
small. Maps from Planck will therefore be dominated by white noise, underneath which there will be low amplitude correlated 
noise associated with striping. Nevertheless, the low amplitude errors from striping introduce a characteristic block-diagonal 
structure in the covariance matrix for the aern- These could confuse searches for small amplitude physical effects, such as 
non-Gaussian features of the CMB signal. 

The implications of our analysis for the application of a hybrid power spectrum estimator to Planck-like data are fairly 
self-evident. Since the noise can be very accurately approximated as white for multipoles i ^ 50, the combination of a number 
of pseudo-Q estimators with different pixel weighting schemes as in E04 should give a close to optimal estimate of the power 
spectrum at high multipoles. At low multipoles, an estimate of Ce can be found by applying a quadratic maximum likelihood 
(QML) estimator (Tegmark 1997c) to a low resolution map. This providesa close to optimal estimate of the power spectrum 
at low multipoles, taking into account of masked regions of the sky, and returns an estimate of the covariance matrix {CiC[) 
that . For the QML estimates, it should be an excellent approximation simply to neglect destriping errors, since these are 
likely to be negligible compared to the cosmic variance. If necessary, destriping errors can be taken into account in the QML 
estimates, and folded into the covariance matrix, by computing the full pixel-noise covariance matrix on a low resolution map 
using an optimal map making algorithm. 

Finally, this analysis has some implications for optimal map making algorithms. The max;imum likelihood map m is given 
by the well-known expression 



where W is orthogonal and A is diagonal. The matrix is a 'prewhitening filter', since the noise matrix of the vector 
n' = W^t is diagonal. In terms of the prewhitened time-stream, equation (27) can be written as. 




(26) 



m = (P'^N"^P)"^P'^N"H, 

where N is the noise covariance matrix Nij = {mnj). Since N is symmetric, it can be written as 
N = WAW'^, 



(27) 



(28) 
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Figure 6. The top panel (a) shows the power spectrum of the input map (solid line) and the power spectrum of the destriped map (filled 
circles). The dashed line shows the power spectrum of the fiducial ACDM model used to generate the input map. The power spectrum 
of the difference (noise) map is also shown. The middle panel (b) shows the power spectrum of the difference map on an expanded scale 
(filled circles) together with the error model shown in Figure 5f. The filled circles in the lower panel (c) show the differences between 
the power spectrum of the destriped map. after subtracting the constant white noise level, and the power spectrum of the input map. 
The solid lines show the dispersion expected from equation (26) using the error model plotted in Figure (6b). The dashed lines show the 
disersion expected from white noise alone. 



P^WA'^W^Pm = P^WA^^W^t. (29) 

Now consider equation (29) applied to a TOD consisting of repeated scans of a single ring with stationary noise of an arbitrary 
power spectral shape. If edge effects are ignored then the components of A will be identical, in which case the solution of 
(29) is obviously the average of t over the rings (equation 4). The rhs of equation (29) effectively filters the TOD, but this is 
undone by the Ihs to return the average of the signal over the rings. 

Now consider a TOD consisting of a set of rings as in the PZancfc-like scanning strategies considered in this paper. The 
maximum likelihood map will differ from a simple average over rings because the rings cross. Thus optimal map making 
performs destriping by comparing the crossing points of the TOD. However, since the noise on a single ring is dominated by 
white noise, optimal map making cannot reduce the errors much below the 'irreducible' white noise errors shown in Figures 
5a and 5d. Optimal map making can, in principle, reduce the map making errors associated with '1//' noise (c/ Figures 5b 
and 5e) if the knee frequency is significantly greater than the spin frequency, but even in this case, the ability to reduce these 
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errors will be limited by the white noise on a ring. These arguments suggest that for a Planck-like scanning strategy'' , simple 
destriping algorithms will be very close to optimal, and may actually be preferable to optimal algorithms because of their 
speed and because thoy require fewer assumptions about the noise. Optimal algorithms arc only optimal if the noise model is 
accurate. In practice, with realistic non-stationary noise, simple destriping may perform just as well and conceivably better 
than an optimal algorithm. 

Acknowledgements: I thank members of CITA for their hospitality during a visit where this work was begun. I am especially 
grateful to Mark Ashdown, Dick Bond, Anthony Challinor and Floor van Leeuwen for helpful discussions. 

REFERENCES 

Bartolo N., Komatsu E., Matarrese S., Riotto A., 2004, submitted to Physics Reports, astro-ph/0406398. 
Bennett, C. et al, 2003, ApJS, 148, 1. 

Brink D.M., Satchlcr G.R,. 1993, Angular Momentum, third edition, Oxford University Press, Oxford. 

Borrill J., Ferreira P. G., Jaffe A. H., Stompor R., 2001, In 'Mining the Sky', Proceedings of the MPA/ESO/MPE Workshop, Edited by 

A. J. Banday, S. Zaroubi, and M.Bartelmann. Springer- Vcrlag, Heidelberg, p403. 
Burigana C, Malaspina M., Mandolesi N., Danese L., Maino D., Bersanelli M., Maltoni M., 1997, Internal Report ITESRE. astro- 

ph/9906360. 

Delabrouille, J., 1998, A&A Suppl. Ser., 127, 555. 

Dore O., Teyssier R., Bouchet F.R., Vibert D., Prunet S., 2001, A&A, 374, 358. 

Efstathiou G., 2004, MNRAS, 439, 603. 

Groth E.J., Peebles P.J.E., 1986, ApJ, 310, 507. 

Hinshaw G. et al, 2003, ApJS, 148, 63. 

Keihanen E., Kurki-Suonio H., Poutanen T., Maino, D., Burigana C., 2003, submitted to A&A. astro-ph/0304411. 
Maino D., et al., 1999, A&A Supph Ser., 140, 383. 

Maddox S.J., Efstathiou G., Sutherland W., 1996, MNRAS, 283, 1227. 
Natoli P., de Gasperis G., Gheller C., Vittorio N., 2001, A&A, 372, 346. 

Poutanen T., Maino, D., Kurki-Suonio H., Keihanen E., Hivon E., 2004, submitted to MNRAS. astro-ph/0404134. 

Revenu B., Kim A., Ansari R., Couchot F., Delabrouille J., Kaplan J., 2000, A&AS, 142, 499. 

Slosar A., Seljak U., 2004, submitted to PRD, astro-ph/0404567. 

Spergel D.N. et al, 2003, ApJS, 148, 175. 

Stompor R., White M., 2004, A&A, 419, 783. 

Tegmark M., 1997a, ApJL, 480, L87. 

Tegmark M., 1997b, PRD, 56, 4514. 

Tegmark M., 1997e, PRD, 55, 5895. 

Tuovinen J., 2003, Planck Newsletter Number 4., p7. (http;//www.rssd.esa.int/SA/PLANCK/does/Newsletters/PlanekNewsletter4.pdf) 
Varshalovich, D.A., Moskalev A.N., Khersonskii V.K., 1988, Quantum Theory of Angular Momentum, World Scientific, Singapore, 
van Leeuwen F., et al, 2002, MNRAS, 331, 975. 
Wright E.L., 1996, astro-ph/9612006. 

Wright E.L., Hinshaw G., Bennett C.L., 1996, ApJ, 458, L53. 



■t The situation is different for a WMAP-type strategy with fast precession (see Bennett et al 2003; Hinshaw et al 2003), in which a 
pixel on the sky is observed on many different timescales. 
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